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Singular gauge fields in the partition function for QCD can lead to a domain-like picture for the QCD vacuum 
by virtue of constraints on quantum fluctuations at the singularities. With a simple model of hyperspherical 
domains with interiors of constant field strength we show that the main features of gluon condensation and an 
area law for static quarks can be realised. The Dirac operator in such a background is exactly soluble. Chirality 
properties of the solutions show agreement with recent lattice results. 


1. INTRODUCTION 

The nature and consequences for confinement 
and chiral-symmetry breaking of long-range gluon 
fields in the vacuum have been the subject of ex¬ 
tended investigation. While many studies rely 
on one or more of the specific configurations of 
monopoles, vortices and instantons, there are ar¬ 
guments 0 that in fact the complete hierarchy 
of singular gauge fields must play a role, espe¬ 
cially if these objects condense in the vacuum. 
It is natural then to seek a model in which the 
bulk properties of this hierarchy can be effectively 
represented. We partially do this here by intro¬ 
ducing a model for such fields in which singu¬ 
larities in vector potentials are concentrated on 
three-dimensional hypersurfaces dVj in Euclidean 
space, in the vicinity of which gauge fields can be 
divided into a sum of a singular pure gauge S^j, 
and a regular fluctuation part and where a 
certain SU (3)coiour colour vector n“ can be asso¬ 
ciated with S^. Thus we have idealised the situ¬ 
ation by ignoring lower dimensional singular ob¬ 
jects which would appear as ’’dislocations” at the 
abelian domain wall singularities. 

The requirement of finiteness of the action 
(density) for singular background plus fluctua¬ 
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tion parts for configurations dominating the func- 
tional integral Q can be formalised in our case as 
specific boundary conditions on dVj for the fluc¬ 
tuation fields charged with respect to rij. The in¬ 
teriors of these regions thus constitute “domains” 
Vj. Gauge field modes neutral with respect to n“ 
are not restricted and generate interactions be¬ 
tween domains. 

Within the present framework then, we decom¬ 
pose a general gauge field = S')) -b Q), and de¬ 
manding finiteness of the classical action density 
at the singularity we come to the conditions 

= 0, a; e dVj, 

for gluons, and 

tjj = = 't[ji jf , X G dVj (1) 

for quarks. The adjoint matrix fij = T°"n'j ap¬ 
pears in the condition for gluons, and a bag-like 
boundary condition arises for quarks, with a unit 
vector if^^x) normal to dVj. 

In a given domain Vj the effect of fluctuations 
in neighbouring regions can be seen as an exter¬ 
nal gauge field neutral with respect to n'j. 
This enables an approximation in which different 
domains are assumed to be decoupled from each 
other but, to compensate, a certain mean field 
is introduced in domain interiors. To make the 
model analytically tractable we consider spherical 
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domains with fixed radius R and approximate the 
mean field in Vj by a configuration with the field 
strength which is (anti-)self- 

dual = ±B\i} so that = B'^Si^p, 

with B a constant the same for all domains. The 
constant colour matrix cos + t® sin 

belongs to the Cartan subalgebra with angles 
€ {f (2fc + 1), A: = 0,..., 5} corresponding to 
the Weyl subgroup. There is no source for this 
field on the boundary and therefore it should be 
treated as strictly homogeneous in all further cal¬ 
culations. The homogeneity itself appears here 
just as an approximation. Thus the partition 
function we will deal with can be written 

1-1 

App 


xexpj-S-^^^^ [qW 


where the gluons and quarks are integrated over 
the spaces Tq and respectively, and the ther¬ 
modynamic limit assumes N oo but with 
the density v~^ = N/V taken fixed and finite. 
The measure of integration over parameters char¬ 
acterising domains is defined as 


/ 


d(7i... 



Here (pi and 6i are spherical angles of the chromo- 
magnetic field, tOi is the angle between the chro- 
momagnetic and chromoelectric fields, is the 
angle in the colour matrix hi, cxi is the chiral an¬ 
gle and Zi is the centre of the domain 14. 

This partition function describes a statistical 
system of density NjV composed of noninter¬ 
acting extended domain-like structures, each of 
which is characterised by a set of internal parame¬ 
ters and whose internal dynamics are represented 
by the fluctuation fields. 


2. MEAN FIELD CORRELATORS 

In this model the connected n—point correlator 
of field strength tensors, 

N 

B%{x) = -{x- zpflR^), 
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can be calculated explicitly using the measure, 
Eq. (|^). Translation-invariant functions 

En{xl,...,Xn) = ^ J d‘^z\\^e{l - ) 

i—1 

emerge and can be seen as the volume of the re¬ 
gion of overlap of n hyperspheres of radius R and 
centres (xi,... ,Xn), normalised to the volume of 
a single hypersphere v = tt^R‘^/2. The functions 
E!„ are continuous and vanish if \xi — Xj\ > 2R. 
Correlations in the background field have finite 
range 2R. The Fourier transform of S„ is then 
an entire analytical function and thus correlations 
do not have a particle interpretation. The statis¬ 
tical ensemble of background fields is not Gaus¬ 
sian since all connected correlators are indepen¬ 
dent of each other and cannot be reduced to the 
two-point correlations. 

The simplest application of the above correla¬ 
tors gives a gluon condensate density which to 
this approximation is = 4i?^. 

Another vacuum parameter which plays a sig¬ 
nificant role in the resolution of the 14^(1) prob¬ 
lem is the topological susceptibility ||] . To define 

this we consider first the topological charge den- 
2 ~ 

sity Q{x) = ■^p 7 Fp,^{x)Fp„{x) which to lowest 
order is 

* 5 ( 2 ;) = COSU3j, 

where luj G {0,7r} depends on the duality of the 
j-th domain. For a given field configuration the 
topological charge is additive Q = fy d'^xQ{x) = 
q(A+ — A_) for —Nq < Q < Nq where q = 
is a ‘unit’ topological charge, namely the 
absolute value of the topological charge of a single 
domain, and A+ (iV_) is the number of domains 
with (anti-)self-dual field, N = A+ -|- A_. The 
probability of finding the topological charge Q in 
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a given configuration is defined by the ratio of the 
number of configurations AfNiQ) with a charge Q 
and the total number of configurations Mn, 

N\ 

^ 2^ {N/2 - Q/2q)l {N/2 + Q/2q)l' 

The distribution is symmetric about Q = 0, 
where it has a maximum for N even. For N odd 
the maximum is at Q = ±q. Averaged topologi¬ 
cal charge is zero. 

The topological susceptibility y is determined 
by the two-point correlator of topological charge 
density, which in the lowest approximation gives 

X = J d'^x{Q{x)Q{0)) = (3) 


3. AREA LAW 


To zeroth order in fluctuations the Wilson loop 
can be written 


N 


W(L) = lim TT 

V J- J- 


(PZi 


V,N^oa Jv V 

.7 = 1 


dai 


N, 


Tiexpii da^y{x)B^v{x) 


'Sl 


where Stokes’ theorem has been used. Note that 
path ordering in our case is not necessary since 
the matrices are assumed to be in the Cartan 
subalgebra. Computationally it is convenient to 
consider a circular contour in the ( 0 : 3 , 374 ) plane 
of radius L with centre at the origin. Details of 
the calculation are given in [^. Put briefly, the 
colour trace is evaluated exactly followed by the 
integration over orientations of the vacuum field. 
Integrating over the domain centres and taking 
the thermodynamic limit {N —>■ 00 , v = V/N = 
7 r^i?^/ 2 ), one obtains for a large Wilson loop L ^ 
R 


W{L) = lim 

N^ao 


1 - 


n N 




with the exponent U{L) = ctttL^ -|- 0{L) display¬ 
ing an area law with string tension for SU (3)coiour 
given by cr = Bf{BR?) with the function 

,, , 2 ( ^ . 


2 V 3 r^^^dx 

-/ — sm X 

TTZ Jo X 

having a purely geometrical origin. This function 
is positive for z > 0 and has a maximum for z = 
1.55. We choose this maximum to estimate the 
model parameters by fitting the string constant 
to the lattice result, 

Vb = 947MeV, R-^ = 760MeV, (4) 

and get for the gluonic parameters of the vac¬ 
uum = 420MeV, X = (197MeV)^ ^{F^} = 
0.081(GeV)"‘ and q = 0.15, while the density of 
the system is 42fm“^. There is no separation of 
characteristic scales y/BR « 1, hence an approx¬ 
imation based on large or small domains is not 
justifiable and has been not used here. 

If B goes to zero then the string constant van¬ 
ishes. This underscores the role of the gluon 
condensate in the confinement of static charges. 
On the other hand, if the number of domains is 
fixed and the thermodynamic limit is defined as 
V,R —> 00 , N = const. < 00 , namely if the do¬ 
mains are macroscopically large, then W{L) = 1, 
which indicates the absence of a linear potential 
between infinitely massive charges in a purely ho¬ 
mogeneous field. 

An area law obviously does not occur for ad¬ 
joint charges because of the presence of zero 
eigenvalues of the adjoint matrix n“T“. 

4. DYNAMICAL CONFINEMENT 

Confinement of the fluctuation fields can be 
seen in the analytical properties of their propaga¬ 
tors. For the above boundary conditions the fluc¬ 
tuation quark and gluon propagators can be an¬ 
alytically calculated by reduction to the charged 
scalar field problem. This in turn is essentially 
just that of a four-dimensional harmonic oscilla¬ 
tor with the orbital momentum coupled to the 
external field, and the general solution has been 
found exactly by decomposition over hyperspher- 
ical harmonics Q]. 

Due to the above boundary conditions the 
X— space propagators of charged fields are defined 
in regions of finite support where they have in- 
tegrable singularities so that their Fourier trans- 
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forms are entire functions in the complex momen¬ 
tum plane. This is consistent with a confinement 
of dynamical fields |^. Entire propagators have 
the physically appealing property of a Regge spec¬ 
trum of relativistic bound states]^. 

5. CHIRALITY OF QUARK MODES 


are known in analytical form and will be dis¬ 
cussed in detail elsewhere [Q . For illustration we 
give the eigenvalue equations for the purely radial 
(fc = 7711 = 7712 = 0) modes. With A = A/\/^, 
2 ; = i?r^/2 we have the equation for A^qqq 

yiAM (A^ -h 1, 3, z) -f 2e^“M (A^, 2, z) = 0 


Turning to quarks specifically, we address the 
eigenvalue problem for the massless Dirac opera¬ 
tor in a domain, 

(ip- A)V'(x) = 0, 

— dpi iBpi — dpi ~fil3pijiXji ^ 

subject to the boundary condition Eq.(|i|). We 
introduce projectors in the colour direction N± = 
i(l ± n/|n|) and spin projection with respect to 
the magnetic field E± = i(l ± 'EB/B) with B = 
\h\B. It is also useful to introduce the mixed 
projector + N-E-c^. Substituting 

tjj = {ip + X)^{x), 

$ = P±$o + RFO-f4>-n + RFO-$-i (5) 
shows that satisfies 
i-D'^ + 2CB + = 0, 

whose exact solution is given in Q. It turns out 
that it is impossible to fulfill the boundary condi¬ 
tion Eq.(|^) with the term $0 in Eq.(||). Thus 4) 
must have definite chirality corresponding to the 
duality of the vacuum gluon field. 

The boundary condition enforces a discretisa¬ 
tion of the eigenvalues A. The states are then 
labelled by the set of quantum numbers: Af = 
77, k , 7771, ?772, s with 77 = 0, 1, ... a radial quan¬ 
tum number labelling the solutions of the bound¬ 
ary condition for given 0(4) = 0(3) x 0(3) an¬ 
gular quantum numbers fc, 7771,7772, C = cor¬ 
responding to the colour projection with respect 
to 77 and s =t,i representing the spin projec¬ 
tion with respect to the magnetic field. Zero 
modes are absent under chosen boundary condi¬ 
tions. The eigenspinors have the structure 

V’TV = ir^X^f + IbXN = TXY, 

where spinors and XN s-s well as the con¬ 
straints defining the corresponding eigenvalues 


and that for A 


i- 

nOOO 


A^ -I- 2 


M (A^ + 3,3,z)-M (A^ -f 2,2, z) 


-ke^*“AM {A^ + 2,2,z) =0 


and AJjqqq — A^qqq and AJ^qqq — A^qqq. Here 
M(a, 5, z) is the confluent hypergeometric func¬ 
tion. The sign minus (plus) in front of ia corre¬ 
sponds to the (anti-)self-dual domain. 

With these solutions we investigate the chiral¬ 
ity of the eigenmodes and also how their probabil¬ 
ity density correlates with the underlying domain. 
Following 1^ we introduce the local chirality pa¬ 
rameter X (x) defined via 


tan(^(l-kA(a:))) 


li’Rix)]' 


which will give extremal values A = ±1 at posi¬ 
tions X where ipix) is purely right(left) handed. 
For lattice overlap fermions, even in uncooled 
gauge fields, histograms of X measured at proba¬ 
bility density maxima for low-lying overlap-Dirac 
eigenmodes show peaks close to the extremal val¬ 
ues 1^. This indicates that low-lying modes are 
strongly chiral and are thus a useful filter for 
the duality of “objects” underlying the gauge 
field fluctuations. Within the present model we 
observe that at the centres of domains all of 
the eigenmodes are exactly chiral and probabil¬ 
ity densities are maximal. The “width” of the 
peaks for the lowest modes at half-maximum 
varies for different values a and is of the order 
of .12 — 0.14fm, consistent with the lattice obser¬ 
vations of As is seen from Fig. the chi¬ 
rality of the lowest mode (77 = 0) monotonically 
decreases with distance from the centre. The chi¬ 
rality parameter for the excited modes alternates 
between extremal values, the number of alterna¬ 
tions is correlated with the radial number 77, and 
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Figure 1. Chirality parameter for the lowest ra¬ 
dial modes 4 ’ii 000 ~'^^ self-dual domain, a = -Kj^. 


the half-width decreases with growing n. The chi¬ 
rality parameter is zero at the boundary for all 
modes. Qualitatively this picture does not de¬ 
pend on angle a. 

The degree of chirality of the modes can be 
conveniently characterized via a histogram if we 
average X(x) over a small neighborhood of the do¬ 
main centre. Then, by assuming that all values 
of a are equally probable we compute the proba¬ 
bility to find a given value of smeared X among 
a set of modes. The result given in Fig.(||) was 
obtained for the lowest modes with all possible 
spin-color orientaions. We observe the double 
peaking close to the extrema with X « ±0.8. 
Including higher modes will broaden the peaks 
and build up a central plateau. This feature as 
well as the above mentioned values for the half¬ 
decay and the density of domains is in qualita¬ 
tive and quantitative agreement with recent lat¬ 
tice results [p|jT^|. However, at this stage, it 
would be premature to go beyond purely nu¬ 
merical comparison of our results and those of 
the lattice. Finally, we quote the estimate of 
the quark condensate density averaged over ^nti- 
)self-dual configurations and a obtained in Q di¬ 
rectly from the quark propagator. For B and R 
as in Eq. (i_ the density at the domain center 
is equal to center = — (228MeV)^. It should 
be stressed that this nonzero result is due to the 
non-integral topological charge q associated with 
domains. This and our results for the chirality 



X 

Figure 2. Histogram of chirality parameter X av¬ 
eraged over the central region with radius 0.05fm. 

parameter X suggest that the underlying mech¬ 
anism of chiral symmetry breaking here is the 
strong chirality of the low-lying non-zero modes, 
which in turn is a consequence of the (anti-)self- 
duality of the background gluonic configurations 
and their non-integral topological charge. 
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